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Stability analysis on the De Sitter universe in pure gravity theory is known to be useful in many 
aspects. We first show how to complete the proof of an earlier argument based on a redundant field 
equation. It is shown further that the stability condition applies to k ^ Friedmann-Robertson- 
Walker spaces based on the non-redundant Friedmann equation derived from a simple effective 
Lagrangian. We show how to derive this expression for the Friedmann equation of pure gravity 
theory. This expression is also generalized to include scalar field interactions. 
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Inflationary theory provides an appealing resolution for the flatness, monopole, and horizon problems of our present 
universe described by the standard big bang cosmology JO . It is known that our universe is homogeneous and isotropic 
to high degree of precision [9. Such a universe can be described by the well known Friedmann- Robertson- Walker 
(FRW) metric There are only three classes of FRW spaces characterized by their topological structure: one can 
either have a closed, open or flat universe according to the observations at large. 

It is also known that gravitational physics should be different from the standard Einstein models near the Planck 
\Q ' scale For example, quantum gravity or string corrections could lead to interesting cosmological consequences 

t-H , [0]. Moreover, some investigations have addressed the possibility of deriving inflation from higher order gravitational 
corrections ||H]. 

A general analysis of the stability condition for a variety of pure higher derivative gravity theories is very useful 
' in many respects. It was shown that a stability condition should hold for any potential candidate of inflationary 
Q\ , universe in the flat FRW space We will first briefly review the approach of reference H based on a redundant 

field equation. The proof will be shown to be incomplete. We will also show how to complete the proof with the help 
q-i of the Bianchi identity for some models where the redundant equation can be recast in a form similar to the Bianchi 
^ • identity in a FRW background. 

In addition, the derivation of the Einstein equations in the presence of higher derivative couplings is known to 
be very complicated. The presence of a scalar field in induced gravity models and dilaton-gravity model makes the 
• ?-h . derivation even more difficult to derive. 

We have developed a simpler derivation by imposing the FRW symmetry before varying the action, while keeping 
the proper time lapse function jl0| . We try to generalize the work in in order to obtain a general formula for the 
non-redundant Friedmann equation. It can be applied to provide an alternative and simplified method to prove the 
validity of the stability conditions in pure gravity theories. In fact, this general formula for the Friedmann equation 
is very useful in many area of interests. 

II. FRIEDMANN EQUATION AND STABILITY OF PURE GRAVITY THEORIES 

The generalized Friedmann-Robertson- Walker (GFRW) metric can be read off directly from the following equation: 
ds 2 = g^ KW dx»dx v = -b(t) 2 dt 2 +a 2 (t)^ T ^^+r 2 dny (1) 
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Here dfl is the solid angle dfl = d9 2 + sin 2 # d\ 2 , and k = 0, ±1 stand for a flat, closed or open universe respectively. 
Note also that the FRW metric can be obtained from the GFRW metric by setting the lapse function bit) equal to 
one, i.e. b = 1, in equation (|l|). 

One can list all non- vanishing components of the curvature tensor as 

R% = \[HB + 2B(H + H 2 )]S), (2) 

RV kl = { H 2 B + ^)C\ v (3) 

Here C l3 kl = e l3m e m ki with e yfc denoting the three space Levi-Civita tensor 0. Here ' denotes differentiation with 
respect to t and H = d/a is the Hubble constant. We have written B = 1/b 2 for later convenience. 

Given a pure gravity model one can cast the action of the system as S = J d A x^fgC = N J dta 3 L(H,H,k/a 2 ) in 
the FRW spaces. Here N is a time independent integration constant. If we take L as an effective Lagrangian, one 
can show that the variation with respect to a gives 

SL , 9 ■ SL d s , , d,SL SL, SL 

3L - H m + - H ^ = {2H+ dt^ H + ^15 + m ] + 2k ~k- (4) 

Note that a 3 L is normally referred to as the effective Lagrangian. We will also call L the effective Lagrangian unless 
confusion occurs. The Above equation is the space-like ij component of the Einstein equation 

G \ iV = t^ v (5) 

with denoting the generalized energy momentum tensor associated with the system. It is known that this equation 
is in fact a redundant equation. Indeed, one can define = G^ v — t^ v and write the field equation as H^ v = 0. 
Hence one has 

DpHV = (6) 

from the energy conservation [D^t^ = 0) and the Bianchi identity {D IJb G tlv = 0). Indeed, the extended Bianchi 
identity (||) can be shown to give 

{dt + 3H)H tt + 3a 2 HH 3 = 0, (7) 

as soon as the FRW metric is substituted into equation (|). Here H 3 = \h l ^Hij and = a 2 hij. It is now 
straightforward to show that i3y = H 3 hij. In fact, equation (7j) indicates that: " H u = implies H 3 = 0" as long as 
a 2 H 0. On the other hand, H 3 — implies instead (dt + 3H)H tt — 0. This implies a 3 H u = constant. Hence the 
ij equation can not imply the Friedmann equation H u = 0. Hence any conclusion derived without the Friedmann 
equation is known to be incomplete. 

We will briefly review the stability analysis obtained from the analysis based on the redundant equation Q 
here and show how to make up the loop-hole in this approach. Suppose that we are given a pure gravity theory, the 
stability of the background inflationary solution for the Hubble constant H = Ha to the redundant field equation (Q) 
can be obtained by perturbing H = Hq + SH. The leading order perturbation equation can be shown to be 

3H F + F = (8) 

along with the zeroth order equation that vanish according to the field equation. This in fact takes a little arguments 
as shown in reference One can show that the zeroth order perturbation equation from the perturbed Friedmann 
equation leads directly to the field equation for the background field. For simplicity the parameter k is set as k = 
in reference ||]. Here F is defined as 

F = L 22 (0)SH + 3H L 22 (0)SH + (6L 2 (0) + 3H L 21 (0) - L n )SH (9) 

In addition, the coefficient of expansions are defined by 

SL SL 

L(H,H) = L(H ,0) + — (H ,0)SH + -^(H o ,0)SH = L(0) +L 1 (0)££f + L 2 {Q)SH (10) 
SH SH 



Y^{H,H) = —(H ,0) + 7TE ^(H ,0)SH + T7777; (H ,0)5H = L 1 (0) + L n (0)5H + L u (0)6H (11) 



SL , ' , SL , s 5 2 L , ^ S 2 L 
W V,B) = m (H ,0) + —^ { H^)6E + 15rH 

SL SL 8^ L S^ L 

j^(H,H) = ^(H ,0) + j^(H ,0)SH+j^(H Q ,0)SH^L 2 (0) + L 2 ^5H + L 22 (0)SH (12) 
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If we focus on the solution F — || , one has 

SH = A+e B+t + A-e B -\ (13) 

Here A± denotes arbitrary constants and 



£ (14) 

denotes the characteristic roots to the characteristic equation 

L 2 2(0)x 2 + 3H L 22 {0)x + 6i 2 (0) + 3H L 21 {0) - L n = (15) 

of the ODE (^). Here A = 9HqL2 2 ~4,L 22 (6L 2 + 3HqL 2 i — Lu) denotes the discriminant of the characteristic equation 
of®. 

One can integrate SH to obtain 

a(t) = a exp ^H t + ^ B+t + l^ 5 "*) ■ ( 16 ) 

Therefore, one finds that stability of the de Sitter type inflationary solution will require both characteristic roots B± 
to be negative. If one of the roots is positive and the other one is negative, then there may exist a limited period 
of inflation. This sort of inflation will come to an end in a time duration of the order of 1/-B p with B p denoting the 
positive root. Choosing a sufficiently small value of l/B p allows inflation to exit naturally |8|]. Therefore the sign 
of the roots to the characteristic equation (|l5|) can be checked to see if the system supports a stable inflationary de 
Sitter solution. If the discriminant is negative, the solution of B± will contain an oscillating phase. Hence the system 
is stable again. Since this argument is based on the redundant field equation, this stability analysis is not complete. 
In other words, the redundant Gij equation will normally take the form of dt{a 3 Gtt) = 0. Hence analysis based on 
the equation will be quite indirect and incomplete. 

There are two problems with this stability condition. First of all, this condition is obtained from the redundant 
equation. One does not know the validity of the field equation, not to mention the stability condition derived 
from it. Secondly, there are homogeneous terms in equation (||) in addition to F = 0, i.e. F = k\ exp[— 3Hot] 
with an arbitrary constant k\. The first problem is not easy to answer for the moment. The second problem can 
be resolved immediately. One notes that the complete solution to the redundant equation (|J) is in fact SH = 
A + e B + t + A-e B - 1 + fci/[(6£ 2 (0) + 3H L 21 (0) - L n )a%\. Here a (t) = cxpH t. This obviously will not affect the 
stability analysis as long as we are interested in the inflationary universe where the particular solution is negligible 
in the above equation unless the denominator of the fti-term happens to vanish. In fact, we are going to show that 
F = is not only a lucky guess, it can be derived from perturbing the Friedmann equation. But one can not be 
sure about this unless a closed form expression for the Friedmann equation is available so that a model independent 
analysis is applicable. 

Nevertheless, one can still resolve this problem by looking into the details of the Bianchi identity. As to the first 
problem with this condition, one notes that in most cases, the redundant equation can be rearranged as 

d t (a 3 H tt ) = (17) 

using the Bianchi identity. The solution to above equation is H tt = constant x a~ 3 . Hence one can show that the 
Friedmann equation has to be of the form 

H tt = F + k ia - 3 = (18) 

if the redundant equation can be written as the combination dt(a 3 F) = with F — the corresponding equation 
leading to the first order equation F = shown in Ref. ||. To be more specifically, 5F = F+ to the leading order 
in SH and its derivatives. Here H = Hq + SH. This follows from the fact that dt[a 3 (H u — F)} = implies that the 
difference H u — F has to be proportional to a~ 3 with some arbitrary constant k\. Therefore, one can effectively work 
with the F = solution if we are working on an inflationary background De Sitter solution. This is because H tt ~ F 
in the De Sitter background. Therefore, any analysis based on the ansatz F = can only be justified in the De Sitter 
background. In particular, stability conditions derived from F = adopted in Ref. H can not be justified from above 
analysis in anti De Sitter space. This is because the undetermined part k\a~ 3 will affect the result significantly. 
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While we suspect that F — should probably be the first order Friedmann equation we are looking for, we are not 
sure if the redundant equation can always be cast into the familiar form shown above. Moreover, the true Friedmann 
equation can look like F + ki /a 3 = even if we can write the redundant equation in above familiar form. Fortunately, 
one can in fact derive a closed form for the Friedmann equation similar to equation (|J). 

The Friedmann equation can be recast as 

d o SL 5L 

L + (H- + 3H>-H) W --H = (19) 

after some algebra. This is done by a variation of £ GFKW with respect to b (or equivalently with respect to gu) and 
setting b = 1 afterwards. Here L GFRW = J d 3 xL(g ilv — <7™ FW ). One notes that the crucial point in the derivation 
is due to the observation that any variation of L with respect to HB has to be equivalent to the variation of L with 
respect to 2BH. This is because the term HB always shows up with 2BH as indicated in the explicit formulae listed 
in equations (§-||). Note that equation ( |l9|) is known as the minimum Hamiltonian constraint TL = nd — L(a, d) = 
in the case where SL/SH — 0. For example, one can write L = —R = 6[k/a 2 — H 2 ] after proper integration by parts. 
Hence the Hamiltonian constraint is identical to the Friedmann equation ( |l9|) in this model. 

Note in particular that even the term containing k does not get involved explicitly in the Friedmann equation, 
equation ( fl9| ) remains valid for arbitrary k. Our derivation leading to equation (|l9|) is based on a pure gravitational 
action. The derivation of the Friedmann equation in the presence of other sources of interactions is straightforward. In 
addition, the Friedmann equation in D-dimensional FRW space ]Tl[ | can also be derived following similar arguments. 

One can then apply the same perturbation, H = Hq+SH , to the Friedmann equation. The zeroth order perturbation 
equation gives exactly the field equation for the background field H = Hq while the leading order in SH gives F = 



identically. Indeed, perturbing equation (19) gives 



L 22 (0)SH + 3H L 22 (0)SH + (6Z 2 (0) + 3H L 21 (0) - L n )SH = (20) 



to the leading order in SH and its derivatives. Note that a-dependent terms always appear in a combination as 
H 2 + k/a 2 in Ft^i as given by equation (|j|). Hence one can ignore the <5a-dependent terms during the inflationary 
phase when H 3> 1/a 2 . This follows from the fact that 5a ~ aSHAt and hence \S{\/ a 2 ) / 5H 2 \ ~ 1/Ha 2 <§; 1 during 
the inflationary phase. Here At is the time duration for the inflationary phase. Therefore, one is able to show that 
the stability conditions in the inflationary phase are indeed given by the result obtained in equation (|l6|). Hence, one 
would never need to worry about any complication that can possibly weaken the validity of the stability condition 
obtained in reference Q . This stability condition hence serves as a screening device for any possible candidates for a 
realistic cosmological universe without any ambiguity. In addition, the stability condition obtained earlier works also 
for curved FRW spaces in the inflationary phase where H ^> 1/a 2 . 

In short, our result states clearly without ambiguity that physically acceptable inflationary models need to meet 
the stability conditions shown earlier in this section. The perturbative stability indicates that a solution with a 
stable mode and an unstable mode can possibly exit the inflationary phase in due time. Our result based on the 
non-redundant Friedmann equation is complete and remains valid for all FRW models. Most of all, working directly 
on the Friedmann equation (19]) we just derived can save us a lot of trouble in the complete analysis. 



III. FRIEDMANN EQUATION FOR SCALAR-GRAVITY THEORY 

The derivation of the Friedmann equation in the presence of a scalar field is in fact rather straightforward since 
complications only arise from complicated curvature interactions. Indeed, the inclusion of scalar interactions intro- 
duces a kinetic term = —^d^dy^g^". It will take the form = ^<fi 2 B(t) if 4>(x,t) — <f){t). Hence the complete 
effective Lagrangian in the GFRW spaces will take the following form 

ba 3 L(g GFRW , 0) = ba 3 L + \a 3 ^ 2 , (21) 

with Lo denoting the graviton Lagrangian plus everything else except the kinetic term of the scalar field X^. Or 
equivalently, Lo = L — with L denoting the complete effective Lagrangian of the theory. Hence, one can show that 
the Friedmann equation becomes 

Lo-T++ (H^ + 3H 2 -H) d -^- S J§H = 0. (22) 
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Note that the minus sign in front of is due to the a 3 6 _1 L combination from the y/g and the g tt component. In 
addition, the variational equation for the (f> field can be directly obtained from the variation of the effective Lagrangian 
L with respect to cj>. 

The method for deriving the Friedmann equation described here can be extended to theories with any form of 
simple gravitational interactions in a straightforward way. For example, one can study the following action with 
Gauss-Bonnet coupling Q 

L = ~R-^d^cj> + f(<P)R% B (23) 
with R 2 GB = R l , l , a i3R lll ' al3 - iR^R^ + R 2 denoting the Gauss-Bonnet term. The effective La grangian is then 

L = 3(H + 2H 2 + A) + I ft + 24(H + H 2 )(H 2 + ^)f(cf>) (24) 
once the FRW metric is applied. The Friedmann equation (^2|) becomes 

3(iJ 2 + A)(l + 8j ff/) = ^ (25) 

Furthermore, the variational equation of <fi is also straightforward. The result is 

df k 
(/) + 3(f>H-24-t-(H + H 2 )(H 2 + ^r)=Q. (26) 
dep a z 

This agrees with the result in ](| while the derivation is much more straightforward. In fact this simple formula for 
the Friedmann equation can also be generalized to any scalar-gravity theory. It helps to reduce the labor in deriving 
gravitational field equations. It is especially helpful when complicated interactions are present and higher derivative 
terms become important. 
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simplified method to prove the validity of the stability conditions in pure gravity theories. This general formula for 
the Friedmann equation is also very useful in many area of interests. 
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